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methylene with N 2 : The role of singlet methylene in prompt NO 


Stephen P. Walch 

Thermosciences Institute, MS 230-3, NASA Ames Research Center, Moffett Field, California 94035 


//j -ya ~ / 'K 


(Received 11 October 1994; accepted 16 June 1995) 




We report calculations of the minimum energy pathways connecting *CH 2 +N 2 to diazomethane and 
diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to 
HCN 2 +H, and of diazomethane to CH 2 N+N. The calculations use complete active space 
self-consistent field (CASSCF) derivative methods to characterize the stationary points and 
internally contracted configuration interaction (ICCI) to determine the energetics. The calculations 
suggest a potential new source of prompt NO from the reaction of ! CH 2 with N 2 to give diazirine, 
and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN 2 , which leads 
to HCN + N( 4 S) on the previously studied CH+N 2 surface. The calculations also predict accurate 0 
K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, 
respectively. © 1995 American Institute of Physics. 
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I. INTRODUCTION 

The reaction of singlet methylene ( 1 CH 2 ) with N 2 along 
with the reaction of CH( 2 I1) with N 2 , which was previously 
studied by Walch, 1 are potentially important in the formation 
of “prompt” NO. 2 While triplet methylene is known to be 
important in some of the reactions leading to prompt NO 
formation, 2 the role of singlet methylene has not been con- 
sidered previously. There have been several previous studies 
of the CH 2 N 2 surface. An early study of the ground state and 
some of the low-lying excited states of diazomethane, at the 
ground state geometry, was carried out by Walch and 
Goddard. 3 More recent work has been carried out by Bold- 
yrev, Schleyer, Higgins, Thomson, and Kramarenko 4 
(BSHTK) who studied the diazirine and diazomethane 
minima and several saddle points on the CH 2 N 2 as well as 
the CHFN 2 and CF 2 N 2 potential energy surfaces. Guimon, 
Khayar, Gracian, Begtrup, and Pfister Guillouzo 5 (GKGBP), 
as part of a study on the decomposition of tetrazole, also 
reported a number of stationary points on the CH 2 N 2 surface, 
including a saddle point connecting diazirine and diaz- 
omethane. These studies were carried out at a lower level of 
theory than that used in the present study and did not provide 
sufficiently detailed information for combustion modelling 
studies. In this paper we report a detailed study of the path- 
ways for addition of l CH 2 to N 2 to give diazomethane and 
diazirine and the pathway for rearrangement of diazirine to 
diazomethane. Calculations with extended basis sets are also 
reported which lead to accurate heats of formation for diaz- 
irine and diazomethane . 

In Sec. II we discuss qualitative features of the reactions 
considered here. Section III contains the technical details of 
the calculations, Sec. IV contains the results, and Sec. V 
concludes the paper. 

II. QUALITATIVE FEATURES 

Figure 1 shows the qualitative features of the orbitals for 
the reaction of ! CH 2 with N 2 via a pathway in which the CH 2 
attacks the tt bond of N 2 . In the l A } ground state of meth- 
ylene ( ] CH 2 ) there is a substantial near degeneracy effect 


between the C 2s lone pair and the empty C 2 p orbital of a” 
symmetry. This leads to a pair of singlet coupled sp hybrid 
orbitals directed above and below the molecular plane, as 
indicated in Fig. 1(a). In the addition of l CH 2 to N 2 , two 
bond pairs undergo major changes. These are the two elec- 
trons in the CH 2 C 2s pair discussed earlier and the in plane 
7r orbital of N 2 . By analogy to the reaction of CH( 2 II) with 
N 2 (Ref. 1) and the reaction of CH( 2 II) with H 2 (Ref. 6), the 
addition of ] CH 2 to N 2 occurs via a pathway, as shown in 
Fig. 1(a), in which the CH and N 2 approach each other in a 



FIG. 1 . The electronic structure of selected stationary points on the potential 
energy surface for CH 2 ( 1 i4 j) +N 2 : (a) sp } ; (b) min 1; (c) 5/? 3 ; (d) distorted 
diazomethane; (e) sp 2 ; (0 niin 2. 
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TABLE I. Computed energies for stationary points on the 'CH 2 + N 2 surface. 


Geometry 

CIl a 

AE b 





(a) c 



Reactants 

- 148.427 21 (-0.469 40) 

0.0 



min l 

-148.471 56(-0.518 97) 

-24.0 



min 2 

- 148.480 05(-0.531 98) 

-33.9 



'ch 2 +n 2 

-148.443 88( -0.484 81) 

-9.7 



Geometry 

CIl a 

AE b 

CI2 a 

AE b 



(b) d 



Reactants 

- 148.391 06(- 0.43026) 

0.0 

- 148.392 69(- 0.431 14) 

0.0 

3 ch 2 +n 2 

148.408 70(-0.44677) 

-10.4 



vdW 

-148.390 96(- 0.43381) 

-0.5 



sp3 

- 148.398 69(-0.44458) 

-4.9 



sp 1 

-148.383 37(-0.42918) 

4.7 

- 148.385 66(-0.430 64) 

4.4 

min 1 

-[48.434 44(— 0.47866) 

-23.3 

-148.435 68(— 0.479 61) 

-23.3 

sp2 



-148.359 65(— 0.405 59) 

19.8 

min 2 

- 148.443 45(- 0.49 145) 

-32.9 


-33.2 

ch 2 n+n 



-148.291 37(-0.332 44) 

64.3 

hcn 2 +h 

-148.282 39(- 0.32 161) 

66.2 



Geometry 

CIl a 

AE b 





<c) C 



Reactants 

- 148.270 49(-0.29892) 

0.0 



min 1 

-148.311 75(-0.34441) 

-21.4 



min 2 

— 148.31934(-0. 35505) 

-29.7 



3 ch 2 +n 2 

- 148.290 72(-0.3 1842) 

-12.3 



Geometry 

CASSCF 

A £ a 





(d) 



Reactants 

-147.998 52 

0.0 



sp l 

-147.967 61 

23.5 



mm 1 

-148.015 41 

-3.5 



spl 

-147.939 63 

40.8 



min 2 

-148.011 55 

-2.7 



vdW 

-148.004 30 

-1.9 



sp3 

-148.001 51 

2.2 




a The energies are in the form ICCI (ICCI + Q+ 148). 

h Relative energies in kcal/mol including zero-point energy (see Table II). 

c cc-pVQZ basis set. 

d cc-pVTZ basis set. 

c cc-pVDZ basis set. 


nearly parallel arrangement. This non-least-motion pathway 
allows the formation of the diazirine molecule [Fig. 1(b)], 
while maintaining large overlaps between the components of 
the two bond pairs which change substantially during the 
course of the reaction (see Ref. 6). By contrast, addition of 
! CH 2 to N 2 via a C 2v pathway requires effectively breaking 
one bond and is a forbidden process. 7 

Another pathway leading to diazomethane involves ini- 
tial formation of a dative bond between one N 2s lone pair of 
N 2 and the empty C 2 p orbital of CH 2 [Fig. 1(c)]. This struc- 
ture evolves to diazomethane [Fig. 1(d)] by simultaneous 
donation of the N 2s pair onto carbon and back donation of 
the C 2s pair into the 7 r system of diazomethane. This pro- 
cess involves no barrier. 

Figure 1(e) shows the orbitals for the saddle point con- 
necting diazirine to diazomethane. The two radical orbitals in 
Fig. 1(e) arise from the CN cr bond of diazirine which is 
being broken. The doubly occupied orbital on the center N 


atom starts as a N 2s -like lone pair, but becomes s/?-like in 
the saddle-point region and eventually becomes a n r orbital 
for diazomethane. The orbitals for diazomethane are shown 
in Fig. 1(f). From Fig. 1(f) it is seen that diazomethane has 
two 7r electrons in a’ orbitals and four 7r electrons in a " 
orbitals (with respect to the molecule plane). The two singly 
occupied a " orbitals shown in Fig. 1(f) are singlet paired as 
a result of through bond coupling from the doubly occupied 
a" orbital on the center N. 3 

III. COMPUTATIONAL DETAILS 

Several different basis sets were used in these calcula- 
tions. For the CASSCF derivative calculations, which were 
used to locate the stationary points, the polarized valence 
double zeta set of Dunning and Hay 8 was used. The basis 
sets for C and N are the ( 9s5p)f[3s2p ] basis augmented by 
a single set of 3d functions with exponent 0.75 for C and 
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TABLE II. Harmonic vibrational frequencies and zero-point energy for ’CH 2 +N 2 . Harmonic frequencies in 
cm -1 obtained from CASSCF calculations with a pVDZ basis set. 


sp 1 vdW mini spl min 2 CH 2 N HCN 2 


Frequencies 
3214 
1951 
1461 
1081 
477 
72 1/ 

3348 

1009 

389 

ZPE a 

0.029 45 
sp3 


3100 

2329 

1509 

534 

97 

75 

3176 

373 

80 

0.025 68 


3307 

1660 

1546 

1033 

1047 

836 

3436 

1166 

1018 

0.034 28 


3313 
1595 
1491 
960 
747 
8 8 8 1 
141 
3458 
1047 

0.029 05 


3360 

1902 

1495 

1162 

560 

328 

3505 

1188 

419 

0.031 71 


3253 

1659 

1438 

1001 

3360 

1021 


0.026 73 


3432 

1570 

612 

928 

1251 

891 


0.019 76 


3142 

2342 

1471 

1180 

258 

248/ 

3240 

1084 

232 


ZPE a 

0.029 50 

a Zero-point energy in E H . The zero-point energy for the reactants is 0.022 97. 


0.80 for N (pVDZ). The H basis is (4s)/{2s] augmented 
with a single set of 2 p functions with exponent 1.00. The 
basis sets used in ICCI calculations are the Dunning correla- 
tion consistent cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets. 9 

The calculations were carried out in C s symmetry with 
the mirror plane in the plane of the paper of Fig. 1. The 
active space for the CASSCF calculations included the eight 
electrons which are depicted in Fig. 1(a). This required six 
active orbitals of a' symmetry and two active orbitals of a" 
symmetry. The resulting CASSCF calculation had eight elec- 
trons distributed among eight orbitals. 

The CASSCF gradient calculations used the SIRIUS/ 
ABACUS system of programs, 10 while the ICCI calculations 
were carried out with MOLPRO. 11,12 All electrons were corre- 
lated except for the C 1 5 and N 1 s like electrons with the 
restriction that no more than two electrons were allowed into 
certain weakly occupied orbitals. These orbitals were the 
highest three a' and highest a orbitals for the calculation 
denoted by CI1 and the highest two a and highest a orbit- 
als for the calculation denoted by CI2. A multireference ana- 
log of the Davidson’s correction 13 was added to the ICCI 
energies and is denoted by + Q . 

IV. DISCUSSION 

The computed ICCI energies obtained at the optimized 
CASSCF geometries are given in Tables 1(a), 1(b), and 1(c) 
for results obtained with the cc-pVQZ, cc-pVTZ, and cc- 
pVDZ basis sets, respectively. The zero-point effects were 
estimated as 1/2 the sum of the harmonic frequencies, which 
are given in Table II, and the zero-point effects are included 
in the relative energies, which are also given in Tables 1(a)— 


1(c). Thus, the relative energies are appropriate for compari- 
son to experimental results corrected to 0 K. The relative 
energies from Table I are also shown in Fig. 2. Table III 
gives the values of selected bond lengths and angles for the 
stationary points obtained here and also gives results from 



FIG. 2. Computed energetics for the reaction of CH 2 ( 1 /4 J ) with N 2 . Note 
that while sp3 is a saddle point at the CASSCF level, it is below the 
reactants energy at the ICCI level and it is probable that there is no barrier 
on this pathway. 
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TABLE III Selected bond distances and angles for stationary points. 


Structure 

r(N-N) 

HN-C) 

r(C-H) 

ZNNC 

ZNCH 

zlHCH 

CAS/pVDZ a 

1.121 

2.743 

vdW 

1.109 

152.4 

86.2 

102.6 

CAS//?VDZ a 

1.158 

1.843 

sp 1 

1.093 

84.2 

105.3 

110.8 

MP2/6-3 1 G* b 

1.256 

min 1 (diazirine) 
1.480 1.083 

64.9 

117.2 

119.4 

SCF/3-21G* c 

1.188 

1.472 

1.070 

65.9 

118.1 


CAS/pVDZ a 

1.253 

1.497 

1.081 

65.3 

117.2 

119.7 

SCF/3-21G* 0 

1.128 

1.408 

sp2 

1.068 

99.5 

116.8 


CAS/pVDZ 3 

1.230 

1.488 

1.081 

109.6 

114.2 

120.9 

MP2/6-3 lG* b 

1.150 

min 

1.311 

2 (di azomethane) 
1.077 

178.6 

116.9 

125.0 

SCF/3-2lG* c 

1.110 

1.293 

1.066 

180.0 

118.4 

123.2 

CAS//?VDZ a 

1.150 

1.320 

1.076 

176.3 

116.4 

125.1 

CAS/p VDZ a 

1.118 

1.850 

sp3 

1.101 

167.2 

93.1 

105.5 

MP2/6-3 lG* b 

1.078 

1.810 

CH 2 -N 2 (/Jj) 

1.091 

163.4 

93.9 

105.4 

MP2/6-3lG* b 

1.111 

1.905 

CH 2 -N 2 (/j 2 ) 

1.076 

66.1 

90.5 

114.8 


“Present work. 
b BSHTK, Ref. 4. 
C GKGBP, Ref. 5. 


the work of BSHTK (Ref. 4) and GKGBP (Ref. 5) for com- 
parison. Table IV gives the Cartesian coordinates for the 
CASSCF optimized structures. 

As discussed in Sec. II, 1 CH 2 plus N 2 can lead to diaz- 
omethane via a pathway involving initial formation of a da- 
tive bond or to diazirine by addition to the N 2 tt bond. The 
pathway leading to diazomethane involves no barrier ( vide 
infra), while the pathway leading to diazirine involves a bar- 
rier of 4 kcal/mol. Diazirine is computed to be about 24 
kcal/mol below J CH 2 plus N 2 . Rearrangement to diaz- 
omethane, which is computed to be about 34 kcal/mol below 
l CH 2 plus N 2 , involves a 20 kcal/mol barrier, with respect to 
reactants (see Fig. 2). 

The stationary point structures for diazirine (min 1), dia- 
zomethane (min 2), and for the van der Waals minimum 
(vDW) between l CH 2 + N 2 , as well as the saddle point struc- 
tures for the addition of l CH 2 to N 2 to give diazomethane 
(sp 3) and diazirine (spl) and the saddle point for rear- 
rangement of diazirine to diazomethane (sp2) are shown in 
Fig. 3. All of these structures are found to be of C s symme- 
try. From Table II it is seen that the minima have all real 
frequencies and the saddle points have only one imaginary 
frequency, which is of a ' symmetry. The HCN 2 structure is 
taken from the work of Ref. 1 . From Table III it is seen that 
in the van der Waals structure the C-N distance is 2.7 A and 
the N 2 is approximately perpendicular to the plane of the 
CH 2 (Z.NCH-86.2). As discussed in Sec. II, this structure 
may be thought of as a datively bonded structure in which 
one of the N 2s lone pairs of the N 2 donates into the empty 


2 p orbital of the CH 2 . Including zero-point energy, the vDW 
structure is bound by 0.5 kcal/mol with respect to 
‘CH 2 +N 2 [see Table 1(b)]. 

Table III also compares the structures obtained here to 
those obtained by BSHTK (Ref. 4) and GKGBP (Ref. 5). 
The structures obtained here for diazirine and diazomethane 
are in good agreement with those obtained by BSHTK at the 
MP2/6-31G* level. One interesting result is that diaz- 
omethane is found to be slightly distorted from C 2v symme- 
try. ziNNC= 176.3 in our calculations compared to 178.6 for 
BSHTK. The agreement with the results of GKGBP is not as 
good, but these calculations are only SCF/3-21G*. BSHTK 
did find a saddle point for conversion of diazomethane to 
diazirine which is qualitatively the same as our sp2 . 

Table III also shows two saddle-point structures 
CH 2 -N 2 (te,) and CH 2 -N 2 (/s 2 ) from BSHTK. 4 CH 2 N 2 (W,) 
is the same as our sp 3 . This saddle-point connects between 
the van der Waals minimum and diazomethane. From Table 
III it is seen that there is good agreement between the 
CASSCF saddle-point structure obtained here and the MP2 
results of BSHTK. CH 2 -N 2 (fs2) is analogous to our sp2 
structure. From Table III it is seen that the geometry obtained 
by BSHTK is qualitatively similar. 

The reaction pathways were characterized by following 
the gradient downhill in both directions away from saddle 
points sp 1 , sp2 , and sp3 . In each case the walk away from 
the saddle point was started by displacing the geometry by 
0.01 times the normal mode corresponding to the negative 
eigenvalue of the Hessian matrix. The geometries along the 


J. Chem. Phys., Vol. 103, No. 12, 22 September 1995 



4934 


Stephen P. Walch: Path for reaction singlet CH 2 with N 2 


TABLE IV. Cartesian coordinates (in a.u.) for the stationary points of 'CH 2 + N 2 . 



X 


Y 

z 

HA 

20.000 000 

Reactants 

1.636 348 

1.317 277 

HB 

20.000 000 


-1.636 348 

1.317 111 

C 

20.000 (XX) 


0.000 000 

0.005 815 

NA 

0.000 000 


0.000 000 

-1.061 015 

NB 

0.000 000 


O.(XX) 000 

1.061 015 

HA 

- 1 .386 720 

vdW 

1.635 625 

-0.254 811 

HB 

-1.386 720 


-1.635 625 

-0.254 811 

CA 

-0.099 095 


0.000 000 

-0.495 841 

NA 

0.152 338 


O.(XX) (XX) 

4.681 562 

NB 

0.792 197 


0.(X)0 (XX) 

6.700 902 

HA 

3.306 430 

sp 1 

1 .700 893 

1 .287 566 

HB 

3.306 430 


-1.700 893 

1 .287 566 

C 

3.156 047 


0.000 000 

0.124 535 

NA 

-0.108 92! 


0.000 000 

-1.088 710 

NB 

-0.659 985 


0.000 000 

1 .029 043 

HA 

3.361 487 

min 1 

1 .765 350 

1 .022 398 

HB 

3.361 487 


- 1 .765 350 

1 .002 398 

C 

2.382 992 


0.000 000 

0.713 766 

NA 

0.288 198 


0.000 000 

-1.188 559 

NB 

-0.424 164 


0.000 000 

1.069 997 

HA 

3.422 894 

spl 

1 .777 934 

0.944 009 

HB 

3.422 894 


- 1 .777 934 

0.944 009 

C 

2.421 472 


0.000 000 

0.837 716 

NA 

0.262080 


0.000 000 

-0.963 517 

NB 

-1.739 341 


0.000 000 

0.217 783 

HA 

3.497 594 

min 2 

1.804 588 

0.792 830 

HB 

3.497 594 


-1.804 588 

0.792 830 

C 

2.560 109 


0.000 000 

0.768 062 

NA 

0.171 030 


0.000 000 

0.052 381 

NB 

-1.946 326 


0.000 000 

-0.436 103 

CA 

-0.051 454 

spl 

0.000 000 

0.055 355 

HA 

-1.309 639 


1.656 983 

-0.006 625 

HB 

-1.309 639 


-1.656 983 

-0.006 625 

NA 

0.090 380 9 


0.000 000 

3.548 880 

NB 

0.640 351 9 


0.000 000 

5.589 015 

HA 

0.000 000 

CH 2 N 

1.774 754 

-1.038 946 

HB 

0.000 000 


-1.774 754 

-1.038 946 

C 

0.000 0(X) 


0.000 000 

0.001 129 

NA 

0.000 000 


0.000 000 

2.434 016 


walk for the pathway connecting ! CH 2 + N 2 to diazirine via 
spl , the pathway connecting diazirine to diazomethane via 
sp 2, and the pathway connecting l CH 2 +N 2 to diaz- 
omethane via sp 3 are shown in Figs. 4, 5, and 6, respec- 
tively. The Cartesian coordinates for the geometries at each 
step along these pathways are given in the appendix (PAPS, 
Ref. 14). 

From Table 1(b) it is seen that sp 3 is below the reactants 
energy, while the CASSCF energy for sp 3 is above the 
CASSCF energy of the reactants. This is interpreted to indi- 
cate that the pathway to diazomethane via sp 3 involves no 
barrier. This zero barrier path to diazomethane is consistent 


with the studies by Milligen and Jacox 15 of the products of 
vacuum-ultraviolet photolysis of methane in an N 2 matrix. 
These authors report the formation of di azomethane even at 
14 K. 

Tables 1(a)- 1(c) give the energetics obtained from the 
ICCITQ calculations using the cc-pVQZ, cc-pVTZ, and the 
cc-pVDZ basis sets, respectively. The binding energies for 
diazirine and diazomethane with respect to ] CH 2 + N 2 were 
estimated by extrapolating the values obtained for the cc- 
pVDZ, cc-pVTZ, and the cc-pVQZ, basis sets to the basis set 
limit using the formula given by Woon. 16 The resulting C-N 
bond dissociation energies are 24.4 kcal/mol for diaz- 
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0 % 9 % 

FIG. 3. Geometries for the stationary points on the l CH 2 + N 2 surface. Up- 
per left sp i ; upper-right min 1; center left sp 2 \ center right min 2; bottom 
left vdW; bottom-right sp ^ . 


omethane and 34.1 kcal/mol for diazirine. These values in- 
clude zero-point effects and are appropriate for 0 K. Using a 
heat of formation of 1 CH 2 of 102.1 kcal/mol obtained from 
the heat of formation of 3 CH 2 of 93 kcal/mol and a singlet- 
triplet splitting of 9.1 kcal/mol gives heats of formation of 
77.7 kcal/mol for diazirine and 68.0 kcal/mol for diaz- 
omethane. Laufer and Okobe 17 report the heat of formation 
of diazirine as 60.6-66 kcal/mol. However, there is an ear- 
lier value of 79.3 kcal/mol by Paulett and Ettinger. 18 The 
discrepancy between our calculated value and the value of 
Laufer and Okobe is surprisingly large and our calculations 
support the value obtained by Paulett and Ettinger. Laufer 
and Okobe also report a heat of formation of 5 1 - 60 kcal/mol 
for diazomethane. 19 This is also surprisingly far from our 
computed value. Calculations for the D e of N 2 using the 
same extrapolation method as used here give an extrapolated 
value of 226.8 kcal/mol compared to an experimental value 
of 228.4 kcal/mol. 20 As a further check on the computed 
results, we also computed the energy of 3 CH 2 +N 2 with each 
of the basis sets. The 3 CH 2 geometry was taken from the 
work of Bauschlicher, Langhoff, and Taylor 21 who carried 
out a definitive study of the singlet-triplet splitting in CH 2 . 
These results are also given in Tables l(a)-I(c). Here it is 


tttt 


FIG. 5. Geometries along the reaction path connecting diazirine with di az- 
omethane. 


seen that the cc-pVQZ basis set gives a separation of 9.7 
kcal/mol compared to 9.1 kcal/mol obtained in Ref. 21 with 
a comparable size ANO basis set. Taking these consider- 
ations together leads us to assign a maximum error of 2 
kcal/mol to our reported 0 K heats of formation for diazirine 
and diazomethane. 

From Ref. 20 it is expected that the barrier heights ob- 
tained with the cc-pVTZ basis set would change by less than 
0.5 kcal/mol if extrapolated to the basis set limit. This occurs 
because the barrier height involves an energy difference be- 
tween a bond that is being broken and a bond that is forming 
and therefore there is a cancellation of errors. Accordingly, 
this extrapolation was not carried out for the barrier heights 
reported in Table 1(b). 

Table 1(d) gives energetics for the stationary points ob- 
tained at the CASSCF level with the pVDZ basis set. As 
expected, these energetics are much less reliable than those 
obtained from the ICCI calculations. The general trends are 
that the binding energies are underestimated and the barrier 
heights are overestimated. In spite of these obvious defects 
in the CASSCF energetics, comparison of CASSCF and 
ICCI stationary points for the N + 0 2 reaction 22 showed that 
the CASSCF method does a very good job of predicting the 
geometry in the directions perpendicular to the reaction co- 
ordinate, but that the saddle-point geometries can be shifted 
somewhat along the direction of the reaction coordinate. An 


Y 


%m 

dm 

dm 

dm 

dm 

dm 

dm 

dm 

dm 


FIG. 4. Geometries along the reaction path connecting 'CH 2 +N 2 with di- 
azirine. 


FIG. 6. Geometries along the reaction path connecting with dia- 

zomethane. 
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extreme example of this effect, in the present calculations, 
occurs for sp 3, where the ICCI method finds an energy at 
the saddle-point geometry below that of the reactants. This 
has been interpreted to mean that the barrier is shifted far out 
into the entrance channel and in fact there is probably no 
barrier. In a case where the barrier persists at the ICCI level, 
the shift along the reaction coordinate is normally much less 
significant. 

We now return to consideration of other product chan- 
nels in the reaction of ’CH 2 with N 2 . Formation of CH 2 N + N 
is endothermic by 64 kcal/mol and involves a singlet-triplet 
crossing. Thus, this pathway is not a likely product channel. 
Dissociation of diazirine to CHN 2 + H can lead to either dou- 
blet or quartet CHN 2 . The quartet state (q min 1 of Ref. 1) is 
17.0 kcal/mol above the doublet state (d min 2). However, 
CHN 2 + H is computed to be endothermic by 66 kcal/mol, 
with respect to reactants, for the doublet state of CHN 2 , and 
is also not a likely direct product. 

More likely processes leading to CHN 2 are hydrogen 
abstraction from diazirine or diazomethane, both of which 
are expected to be important products of the reaction of 
'CH 2 with N 2 . Hydrogen abstraction from diazomethane 
leads to an HCN 2 geometry like the dative structure of 
Manaa and Yarkony. 23 Previous theoretical studies 1 indicate 
that there is a large barrier between the dative structure and 
the region of the surface leading to HCN + N. Thus, this path- 
way is not expected to contribute to prompt NO formation. 
On the other hand, hydrogen abstraction from diazirine leads 
to an HCN 2 geometry like that involved in the doublet- 
quartet surface crossing region of the CH + N 2 surface. For 
example, the reaction 

diazirine + OH— >CHN 2 + H 2 0 (1) 

is exothermic by 28.5 kcal/mol. In addition, assuming no 
vibrational relaxation, diazirine formed from 1 CH 2 -fN 2 
would have about 38 kcal/mol of internal energy (energy of 
sp l minus energy of min 2), leading to a total of 66 kcal/mol 
of available energy. As the NN bond is considerably shorter 
in diazirine it is probable that CHN 2 will be formed with a 
large amount of energy in the NN stretching mode. [The 
quartet state of CHN 2 has a longer NN bond (2.377 A) than 
the doublet state (1.635 A), which in turn is longer than the 
NN bond in diazirine (1.253 A).] From Ref. 1 the saddle 
points leading to CH + N 2 on the doublet surface and to 
HCN + N on the quartet surface are 31 kcal/mol and 36 kcal/ 
mol above d min 2, respectively. Thus, if this much energy is 
transfered into vibrational energy, these channels are ener- 
getically accessible. Furthermore, energy in the NN stretch- 
ing mode also favors the doublet to quartet crossing [In Eq. 
(1) the CHN2 is formed initially in the doublet state], leading 
to HCN + N( 4 S), as for the CH( 2 I1) + N 2 reaction. 

V. CONCLUSIONS 

The reaction pathways for the reaction of *CH 2 with N 2 
have been studied using complete active space self- 
consistent field (CASSCF) derivative methods to character- 
ize the stationary points and internally contacted configura- 
tion interaction (ICCI) to determine the energetics. 


] CH 2 +N 2 can lead to diazomethane, with no barrier, or to 
diazirine, with a barrier of 5 kcal/mol. Diazirine is computed 
to be about 24 kcal/mol below ’CH;, plus N 2 . Rearrangement 
to diazomethane involves a 20 kcal/mol barrier, with respect 
to reactants, and diazomethane is computed to be about 34 
kcal/mol below l CH 2 plus N 2 . The CH 2 N+N and CHN 2 +H 
channels are computed to be endothermic by 64 kcal/mol and 
66 kcal/mol, respectively. Thus, these product channels are 
unlikely to be important. 

The computed potential energy surface would suggest 
that diazirine will be readily formed from ] CH 2 plus N 2 . It is 
proposed that subsequent reaction of diazirine with hydrogen 
abstracters such as OH, which would be abundant in the 
combustion environment, would lead to doublet CHN 2 , with 
considerable vibrational excitation, especially in the NN 
stretch mode, which would facilitate doublet to quartet cross- 
ing, leading to HCN + N product. Thus, this pathway may be 
an important new source of “prompt 1 ’ NO. 
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